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Figure 5. Quenched and dynamical decay con- 
stants vs. lattice spacing. 



4. CONCLUSIONS 

To really confront the question of whether sea 
quarks are important we need to push to smaller 
values of the sea quark mass. We also need either 
to push to smaller values of the lattice spacing 
or to continue to develop techniques which allow 
one to carry out simulations at large lattice spac- 
ing which have smaller intrinsic discretization sys- 
tematics than present simulations do. 

This work was supported by the U.S. De- 
partment of Energy and by the National Sci- 
ence Foundation. Simulations were performed at 
the Supercomputer Computations Research Insti- 
tute. 
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Bernard, Labrenz, and Soni [9], Kronfeld [10], 
and Mackenzie [11] argue that the appropriate 
quark mass at which the matrix element is mea- 
sured is not mi but rri2 since it enters in the ki- 
netic energy while m\ is just an overall additive 
constant. Their analysis suggests that we correct 
for this error by adjusting the meson mass. 



aM — ► aM' = aM + (arri2 — ami). 



(6) 



This is a shift of no more than 0.125 at k = 
0.1390. 

Fig. 3 displays a plot of afp\/aMp vs 1/aMp 
for heavy-light systems, including the extrapo- 
lated zero light quark mass points. These are 
k = 0.1675 simulations where crosses, diamonds, 
and squares show the local operators and oc- 
tagons, bursts and fancy diamonds show the non- 
local operators. The fancy squares show the ex- 
trapolation to K c . 

Finally, we convert to physical units by fixing 
the lattice spacing from f w . We plot fp\/Mp vs 
l/Mp from our dynamical Wilson data in Fig. 
4. Diamonds and octagons are local and nonlocal 
currents at k = 0.1670; squares and crosses are 
the corresponding currents at k = 0.1675. 

We calculate f D = 215 ± 40 ± 35 ± 5 MeV 
and (from a long extrapolation to the B mass) 
j B = 150 ± 40 ± 40 ± 5 MeV. The three uncer- 
tainties represent lattice spacing, choice of oper- 
ator (the local operators give about 250 MeV for 
jo and the nonlocal ones, about 180 MeV), and 
a s uncertainty. The main sources of error are 
all systematic. Statistical errors for a calculation 
with a particular operator and choice of pertur- 
bative correction are never more than 5-10 MeV. 
We also find j Ds = 288 ± 45 ± 45 ± 5 MeV (The 
local operator gives 330 MeV, the nonlocal oper- 
ator, 246 MeV.) We have included no uncertainty 
associated with the sea quark mass; it is lumped 
in with the statistical/extrapolation uncertainty. 
We do not see any observable effects of different 
sea quark masses. These calculations are a bit 
higher than quenched calculations done at smaller 
values of the lattice spacing [9,12]. 

It may be that the effect of dynamical fermions 
is to push up the matrix elements but it is also 
possible (and more likely, in our opinion) that 
the large lattice spacing induces a systematic shift 



upwards in the decay constant, especially for the 
tadpole-improved matrix elements. Fig. 5 shows 
a comparison of decay constants from our work 
and from the quenched simulations of Ref. [9], 
using the local axial current in both cases. The 
reader is invited to draw his/her/its own conclu- 
sions! 

0.3 i 1 1 1 1 1 1 1 1 1 1 ,— 
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Figure 3. Lattice pseudoscalar decay constants. 
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Figure 4. Pseudoscalar decay constants extrapo- 
lated to the continuum. 
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plaquette 

-\n(^TrUp) = 

4.18879a s (3.41/a)[l- (1.185 ■ 
0(af)] 



0.070ri/)a s 



(3) 



and is run down to a scale 0(l/a). These are 
our conventions. We have already published the 
uncorrected lattice data (as well as non-tadpole- 
improved results) so that if the theoretical sit- 
uation changes one can go back and recompute 
everything. 

Finally, if we are not at zero mass, 
/(m) = ^(1 - 0.75ki/k c )(1 - 0.75k 2 /k c )- This 
is the "exp(ma)" field renormalization factor of 
Bernard, Labrenz and Soni [9]. 
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Figure 1. Sea quark mass dependence of vector 
meson decay constant. 



3. SOME RESULTS 

3.1. Vector currents 

Fig. 1 shows 1/ fv the vector meson decay con- 
stant using the conserved (Wilson) current from 
our dynamical Wilson simulations (squares with 
k = 0.1670, diamonds with k = 0.1675). Fig. 
2 shows three different lattice currents (crosses 
for the Wilson current, diamonds for the point- 
split current and squares for the local current) 
all at k = 0.1675, all with tadpole improvement. 
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Figure 2. Operator dependence of vector meson 
decay constants. 



The effects of sea quarks are much smaller than 
the systematic differences between different lat- 
tice choices for the operator. 

3.2. Axial currents 

Our extraction of pseudoscalar decay constants 
parallels other recent quenched analyses of these 
quantities. Note however that our lattice spac- 
ing is considerably larger than what is used in 
contemporary quenched simulations. This intro- 
duces an unknown systematic effect on our re- 
sults. 

In order to arrive at physical numbers we then 
carried out the following steps: First, we extrap- 
olated heavy meson decay constants to zero light 
quark mass by a linear extrapolation in the light 
quark hopping parameter to k c , using the two 
lightest quark hopping parameters in each data 
set (k = 0.1615 and either 0.1670 or 0.1675). This 
extrapolation included the K-dependent field nor- 
malization and appropriate Z-factor. Our heavy 
quarks have masses which are large compared to 
an inverse lattice spacing. In this limit the dis- 
persion relation for free Wilson fermions is 



E(k) = m 1 + + . . . 

2m 2 

where m\ = log(/« c /(2/«) — 3) and 
exp(ami) sinh(ami) 



arri'z 



1 + sinh(ami) 



(4) 



(5) 



1 



Simple Matrix Elements with Dynamical Fermions 

The High Energy Monte Carlo Grand Challenge[ *] 

We report on studies of simple matrix elements from simulations with two flavors of sea quarks, both staggered 
and Wilson. We show the decay constants of vector and pseudoscalar mesons. The effects of sea quarks are small. 
These simulations are done at relatively large lattice spacing compared to most quenched studies. 



1. INTRODUCTION 

The HEMCGC has recently completed a se- 
ries of simulations of QCD with two flavors of 
dynamical fermions on 16 3 x 32 lattices at a 
gauge coupling [3 = 5.6 with dynamical staggered 
fermions[l] and dynamical Wilson fermions[2] at 
[3 = 5.3. Here we would like to focus on our cal- 
culations of matrix elements and the effects of sea 
quarks on them. Some of this work has also been 
published in Ref. [3]. 

We wish to update one aspect of the Wilson 
spectroscopy: at Lattice 92 we[4] reported dis- 
crepancies in hadron masses calculated using in- 
terpolating fields on different source timeslices of 
the lattice. These effects have largely gone away 
with more statistics; we believe that they were 
due to long simulation time correlations in the 
data which led to an underestimate of uncertain- 
ties. 

In his summary talk Weingarten[5] remarked 
that our time correlations were so long that we 
probably only had ten uncorrelated lattices in 
our data set. This is a gross misrepresentation 
of our data. Our autocorrelations are so long we 
probably only have two uncorrelated lattices. We 
encourage anyone who is going to revisit Wilson 
spectroscopy at [3 = 5.3 to collect three to five 
times as much data as we did, to be certain that 
time correlations in the data are under control. 

2. OUR VERSION OF TADPOLE IM- 
PROVEMENT (NOVEMBER 1993) 

The lattice is a UV regulator and changing 
from the lattice cutoff to a continuum regulator 
(like MS) introduces a shift 

</|0-* (/U = l/ a )|i)— = 



a D (l + A a.) + . . :)(f\O latt (a)\i)f(m) 

+0(a) + ... (1) 

where A = ^r(Cj^- Q a tt), f(m) converts the 
lattice field renormalization to the continuum and 
the factor aP converts the dimensionless lattice 
number to its continuum result. The 0(a) cor- 
rections arise because the lattice operator is not 
"improved". The C's are calculable in perturba- 
tion theory. Lepage and Mackenzie[6] have pro- 
posed a method for achieving a more convergent 
perturbation expansion called tadpole improve- 
ment, which we have attempted to implement in 
our calculations. Unfortunately there does not 
appear to be universal agreement on how to do 
this for all operators. Our present implementa- 
tion of tadpole improvement (not the one we have 
used in Refs. [2] and [3]) is as follows: 

For massless quarks and local operators Aq 
consists of a sum of two terms: Aq = Apt + A K 
where Apt are the one loop perturbative correc- 
tions computed by (for example) Groot, Hoek and 
Smit[7] or Martinelli and Zhang[8]. A K absorbs 
the perturbative shift in k c , A k = 1.3643 

J- = ii l-A K a s ) (2) 

and is included while changing the "conventional" 
bilinear field renormalization from 2k c to 1/4. For 
"point-split" operators ipTU ip tadpole improve- 
ment makes the substitution U —> uo(U/uo) and 
one is supposed to expand perturbatively in U /«o- 
The quantity u = (\TrU P fl A = 1 - 1.0472a s 
and so there is an additional contribution to 
A = Apt + 1.3643 - 1.0472 for nonlocal op- 
erators regardless of spatial orientation. 

The coupling constant is defined through the 



